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BACKGROUND OF THE INVENTION 
The past years have seen a dynamic change in the ability of science to comprehend 
vast amounts of data. Pioneering technologies such as nucleic acid arrays allow scientists 
to delve mto the world of genetics in far greater detail than ever before. Exploration of 
genomic DNA has long been a dream of the scientific community. Held within the 
complex structures of genomic DNA lies the potential to identify, diagnose, or treat 
diseases like cancer, alzheimers or alcoholism. Answers to the world's food distribution 
problems may be held within the exploitation of genomic information fi-om plants and 
animals. 

It is estimated that by the Spring of 2000 a reference sequence of the entire human 
genome will be sequenced allowing for types of genetic analysis that were never before 
possible. Novel methods of sample preparation and sample analysis are needed to provide 
for the fast and cost effective exploration of complex samples of nucleic acids, particularly 
genomic DNA. 

SUMMARY OF THE INVENTION 
The present invention provides a flexible and scalable metiiod for analyzing 
complex samples of nucleic acids, such as genomic DNA. These methods are not limited 
to any particular type of nucleic acid sample: plant, bacterial, animal (mcluding human) 
total genome DNA, RNA, cDNA and the like may be analyzed using some or all of the 
methods disclosed in this invention. The word "DNA" may be used below as an example 


of a nucleic acid. It is understood that this term includes all nucleic acids, such as DNA 
and RNA, unless a use below requires a specific type of nucleic acid. This invention 
provides a powerful tool for analysis of complex nucleic acid samples. From 
experimental design to isolation of desired fragments and hybridization to an appropriate 
array^ the invention provides for faster, more efficient and less expensive methods of 
complex nucleic acid analysis. 

The present invention provides for novel methods of sample preparation and 
analysis comprising managing or reducing, in a reproducible maimer, the complexity of a 
nucleic acid sample. The present invention eliminates the need for multiplex PGR, a time 
intensive and expensive step in most large scale analysis protocols, and for many of the 
embodiments the step of complexity reduction may be performed entirely in a single tube. 
The invention further provides for analysis of the sample by hybridization to an array 
which may be specifically designed to interrogate jfragments for particular characteristics, 
such as, for example, the presence or absence of a polymorphism. The invention further 
provides for novel methods of using a computer system to model enzymatic reactions in 
order to determine experimental conditions and/or to design arrays. In a preferred 
embodiment the invention discloses novel methods of genome -wide polymorphism 
discovery and genotyping. 

In one embodiment of the invention, the step of complexity management of the 
nucleic acid sample comprises enzymatically cutting the nucleic sample into fragments, 
separating the fragments and selecting a particular fragment pool. Optionally, the selected 
fragments are then ligated to adaptor sequences containing PGR primer templates. 

In a preferred embodiment, the step of complexity management is performed 
entirely in a single tube. 

In one embodiment of complexity management, a type lis endonuclease is used to 
digest the nucleic acid sample and the fragments are selectively ligated to adaptor 
sequences and then ampUfied. 


In another embodiment, the method of complexity management utilizes two 
restriction enzymes with different cutting sites and frequencies and two different adaptor 
sequences. 

In another embodiment of the invention, the step of complexity management 
comprises performing the Arbitrarily Primed Polymerase Chain Reaction (AP PGR) upon 
the sample. 

In another embodiment of the invention, the step of complexity management 
comprises removing repeated sequences by denaturing and reannealing the DNA and then 
removing double stranded duplexes. 

In another embodiment of the invention, the step of complexity management 
comprises hybridizing the DNA sample to a magnetic bead which is bound to an 
oligonucleotide probe containing a desired sequence. This embodiment may further 
comprise exposing the hybridized sample to a single strand DNA nuclease to remove the 
single stranded DNA, Ugating an adaptor sequence containing a Class II S restriction 
enzyme site to the resulting duplexed DNA and digesting the duplex with the appropriate 
Class II S restriction enzyme to release the magnetic bead. This embodiment may or may 
not comprise amplification of the isolated DNA sequence. Furthermore, the adaptor 
sequence may or may not be used as a template for the PCR primer. In this embodiment, 
the adaptor sequence may or may not contain a SNP identification sequence or tag. 

In another embodiment, the method of complexity management comprises 
exposing the DNA sample to a mismatch binding protein and digesting the sample with a 
3* to 5' exonuclease and then a single strand DNA nuclease. This embodiment may or may 
not include the use of a magnetic bead attached to the mismatch binding protein. 

BRIEF DESCRIPTION OF THE FIGURES 
Figure 1 is a schematic representation of a method of complexity management 
comprising restriction enzyme digest, fragment separation, and isolation and purification 
of a fragment size range of interest. 
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Figure 2 is a schematic representation of a method of complexity management 
comprising restriction enzyme digest, fragment separation, isolation and piorification of a 
fragment size range of interest, ligation of an adaptor sequence to the desired fragments 
and amplification of those fragments. 
5 Figure 3 depicts the effect on complexity of PGR amplification using primers with 

and without specific nucleotides. 

Figure 4 is a schematic representation of a method of complexity management 
comprising a type lis restriction enzyme digest, adaptor sequence ligation and 
amplification of desired fragments. 

10_^ 

j Figure 5 depicts type lis restriction enzymes and their cleavage sites. 

;2 Figure 6 is a schematic representation of a method of complexity management 

.C comprising a type lis restriction enzyme digest, adaptor sequence ligation and 
5 !^ amplification of desired fragments. 
15D Figure 7 is a schematic representation of a method of complexity management 

O comprising AP PGR. 

^2 Figure 8 depicts the results of AP PGR on human genomic DNA. 

lU Figure 9 depicts the reproducibility of AP PGR. 

J T Figxare 10 is a schematic representation of a method of complexity management 

20 comprising removing repetitive sequences by denaturing and reannealing genomic DNA. 

Figure 1 1 is a schematic representation of a method of complexity management 
comprising hybridizing a probe sequence attached to a magnetic bead to a pool of 
fractionated DNA. 

Figure 12 is a schematic representation of a method of complexity management 
25 comprising hybridizing a probe sequence bound to a magnetic bead to a pool of 

fractionated DNA, ligating an adaptor sequence containing a class lis restriction enzyme 
site to the DNA/probe duplex, digesting the duplex, ligating a second adaptor sequence to 
the duplex and amplifying. 


Figure 13 is a schematic representation of a method of complexity management 
comprising hybridizing a probe sequence bound to a magnetic bead to a pool of 
fractionated DNA, ligating an adaptor sequence containing a class lis restriction enzyme 
site to the DNA/probe duplex, digesting the duplex, ligating a second adaptor sequence to 
the duplex and amplifying. 

Figure 14 depicts a chimeric probe array. 

Figure 15 is a schematic representation of a method of complexity management 
comprising hybridizing a probe sequence attached to a magnetic bead to a pool of 
fractionated DNA, ligating an adaptor sequence containing a class lis restriction enzyme 
site to the DNA/probe duplex, digesting the duplex, ligating a second adaptor sequence to 
the duplex, amplifying and hybridizing the amplicons to a chimeric probe array. 

Figure 16 is a schematic representation of a method of complexity management 
comprising hybridizing a mismatch binding protein to DNA containing a polymorphism 
and isolating the region containing the polymorphism. 

Figure 17 is a schematic representation of a method of complexity management 
comprising attaching a magnetic bead to the mismatch binding protein of Figure 16. 

Exhibit 1 is an example of one type of computer program which can be written to 
model restriction enzyme digestions. 

Exhibit 2 is an example of one type of computer program which can be written to 
model ligation reactions. 

DETAILED DESCRIPTION OF THE PRESENT INVENTION 
This application relies on the disclosure of other patent applications and literature 
references. These documents are hereby incorporated by reference in their entireties for all 
purposes. 

Definitions 

A "genome" is all the genetic material in the chromosomes of an organism. DNA 
derived from the genetic material in the chromosomes of a particular organism is genomic 



DNA. A genomic library is a collection of clones made from a set of randomly generated 
overlapping DNA fragments representing the entire genome of an organism. 

An "oligonucleotide" can be nucleic acid, such as DNA or RNA, and single- or 
double-stranded. Oligonucleotides can be naturally occurring or synthetic, but are 
5 typically prepared by synthetic means. Oligonucleotides can be of any length but are 

usually at least 5, 10, or 20 bases long and may be up to 20, 50, 100, 1,000, or 5,000 bases 
long. A polymorphic site can occur within any position of the oligonucleotide. 
Oligonucleotides can include peptide nucleic acids (PNAs) or analog nucleic acids. See 
US Patent Application No. 08/630,427 filed 4/3/96. 
10 An array comprises a solid support with nucleic acid probes attached to said 

,Q support. Arrays typically comprise a plurality of different oligonucleotide probes that are 
;2 coupled to a surface of a substrate in different known locations. These arrays, also 

described as "microarrays" or colloquially "chips" have been generally described in the art, 
for example, U.S. Pat. Nos. 5,143,854, 5445934, 5,744,305, 5,677,195 and PCX Patent 
ifi Publication Nos. WO 90/1 5070 and 92/1 0092. Each of which is incorporated by reference 
13 in its entirety for all purposes. These arrays may generally be produced using mechanical 
. "I synthesis methods or light directed synthesis methods which incorporate a combination of 
! y photolithographic methods and solid phase synthesis methods. See Fodor et al., Science, 
il 251:161-111 (1991), Pirrung et al, U.S. Pat. No. 5,143,854 (see also PCT Application No. 
20 WO 90/15070) and Fodor et al., PCT Publication No. WO 92/10092 and U.S. Pat. No. 
5,424,186, each of which is hereby incorporated in its entirety by reference for all 
purposes. Techniques for the synthesis of these arrays using mechanical synthesis methods 
are described in, e.g., U.S. Pat. No. 5,384,261, incorporated herein by reference in its 
entirety for all purposes. Although a planar array surface is preferred, the array may be 
25 fabricated on a STorface of virtually any shape or even a multiplicity of surfaces. Arrays 

may be nucleic acids on beads, fibers such as fiber optics, glass or any other appropriate 
substrate, see US Patent Nos. 5,770,358, 5,789,162, 5,708,153 and 5,800,992 which are 
hereby incorporated in their entirety for all purposes. Arrays may be packaged in such a 
maimer as to allow for diagnostics or other manipulation of in an all inclusive device, see 
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for example, US Patent Nos. 5,856,174 and 5,922,591 incorporated in their entirety by 
reference for all purposes. 

Hybridization probes are oligonucleotides capable of binding in a base-specific 
manner to a complementary strand of nucleic acid. Such probes include peptide nucleic 
5 acids, as described in Nielsen et al, Science 254, 1497-1500 (1991), and other nucleic acid 
analogs and nucleic acid mimetics. See US Patent Application No. 08/630,427 filed 
4/3/96. 

Hybridizations are usually performed under stringent conditions, for example, at a 
salt concentration of no more than 1 M and a temperature of at least 25 °C. For example, 
10 conditions of 5X SSPE (750 mM NaCl, 50 mM NaPhosphate, 5 mM EDTA, pH 7.4) and a 

■ 2 temperature of 25-30°C are suitable for alleie-specific probe hybridizations. For stringent 
^0 conditions, see, for example, Sambrook, Fritsche and Maniatis. "Molecular Cloning A 
,C laboratory Manual" 2"^ Ed. Cold Spring Harbor Press (1989) which is hereby incorporated 
"r. by reference in its entirety for all purposes above. 

Bil Polymorphism refers to the occurrence of two or more genetically determined 

altemative sequences or alleles in a population. A polymorphic marker or site is the locus 
at which divergence occurs. Preferred markers have at least two alleles, each occurring at 
i U fi-equency of greater than 1%, and more preferably greater than 10% or 20% of a selected 
J I population. A polymorphism may comprise one or more base changes, an insertion, a 
20 repeat, or a deletion, A polymorphic locus may be as small as one base pair. Polymorphic 
markers include restriction fragment length polymorphisms, variable number of tandem 
repeats (VNTR's), hypervariable regions, minisatellites, dinucleotide repeats, trinucleotide 
repeats, tetranucleotide repeats, simple sequence repeats, and insertion elements such as 
Alu. The first identified allelic form is arbitrarily designated as the reference form and 
25 other allelic forms are designated as altemative or variant alleles. The allelic form 

occurring most frequently in a selected population is sometimes referred to as the wildtype 
form. Diploid organisms may be homozygous or heterozygous for allelic forms, A 
diallelic polymorphism has two forms. A triallelic polymorphism has three forms. 



A single nucleotide polymorphism (SNP) occurs at a polymorphic site occupied by 
a single nucleotide, which is the site of variation between allelic sequences. The site is 
usually preceded by and followed by highly conserved sequences of the allele (e.g., 
sequences that vary in less than 1/100 or 1/1000 members of the populations). 

A single nucleotide polymorphism usually arises due to substitution of one 
nucleotide for another at the polymorphic site. A transition is the replacement of one 
purine by another purine or one pyrimidine by another pyrimidine. A transversion is the 
replacement of a purine by a pyrimidine or vice versa. Single nucleotide polymorphisms 
can also arise from a deletion of a nucleotide or an insertion of a nucleotide relative to a 
reference allele. 

An individual is not limited to a human being, but may also include other 
organisms including but not limited to mammals, plants, bacteria or cells derived from any 
of the above. 

General 

The present invention provides for novel methods of sample preparation and 
analysis involving managing or reducing the complexity of a nucleic acid sample, such as 
genomic DNA, in a reproducible manner. The invention further provides for analysis of 
the above sample by hybridization to an array which may be specifically designed to 
interrogate the desired fragments for particular characteristics, such as, for example, the 
presence or absence of a polymorphism. The invention further provides for novel methods 
of using a computer system to model enzymatic reactions in order to determine 
experimental conditions before conducting any actual experiments. As an example, the 
present techniques are useftil to identify new polymorphisms and to genotype individuals 
after polymorphisms have been identified. 

Generally, the steps of the present invention involve reducing the complexity of a 
nucleic acid sample using the disclosed techniques alone or in combination. None of these 
techniques require multiplex PGR and most of them can be performed in a single tube. 
With one exception (AP PGR), the methods for complexity reduction involve fragmenting 
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the nucleic acid sample, often, but not always by restriction enzyme digest. The resulting 
fragments, or in the case of AP PGR, PGR products, of interest are then isolated. The 
isolation steps of the present invention vary but may involve size selection or direct 
amplification, often adaptor sequences are employed to facilitate isolation. In a preferred 
5 embodiment the isolated sequences are then exposed to an array v^hich may or may not 

have been specifically designed and manufactured to interrogate the isolated sequences. 
Design of both the complexity management steps and the arrays may be aided by the 
computer modeling techniques v^hich are also described in the present invention. 

10^ Complexity management 

The present invention provides for a number of novel methods of complexity 
management of nucleic acid samples such as genomic DNA. These methods are disclosed 
below. 

1 2 A number of methods disclosed herein require the use of restriction enzymes to 

15J fragment the nucleic acid sample. Methods of using a restriction enzyme or enzymes to 
□ cut nucleic acids at a large rmmber of sites and selecting a size range of restriction 
,2 fragments for assay have been shown. This scheme is illustrated in Figure 1. 
W In one embodiment of the im^ention, schematically illustrated in Figure 2, 

u restriction enzymes are used to cut the nucleic acids in the sample (Fig. 2, Step 1). In 
20 general, a restriction enzyme recognizes a specific nucleotide sequence of four to eight 

nucleotides (though this number can vary) and cuts a DNA molecule at a specific site. For 
example, the restriction enzyme Eco RI recognizes the sequence GAATTC and will cut a 
DNA molecule between the G and the first A. Many different restriction enzymes are 
known and appropriate restriction enzymes can be chosen for a desired result. For 
25 example, restriction enzymes can be purchased from suppliers such as New England 

Biolabs. Methods for conducting restriction digests will be known to those of skill in the 
art, but directions for each restriction enzyme are generally supplied with the restriction 
enzymes themselves. For a thorough explanation of the use of restriction enzymes, see for 


example, section 5, specifically pages 5.2-5.32 of Sambrook, et aL, incorporated by 
reference above. 

After restriction enzyme digestion, the method further requires that the pool of 
digested DNA fragments be separated by size and that DNA fragments of the desired size 

5 be selected (Figure 2, step 2) and isolated (Figure 2, Step 3). Methods for separating DNA 
fragments after a restriction digest will be well known to those of skill in the art. As a non- 
limiting example, DNA fragments which have been digested with a restriction enzyme 
may be separated using gel electrophoresis, see for example, Maniatis, section 6. In this 
technique, DNA fragments are placed in a gel matrix. An electric field is applied across 

10 _ the gel and the DNA fragments migrate towards the positive end. The larger the DNA 

fragment, the more the fargment's migration is inhibited by the gel matrix. This allows for 
the separation of the DNA fragments by size. A size marker is run on the gel 
simultaneously with the DNA fragments so that the fragments of the desired size may be 
hi identified and isolated from the gel. Methods for purification of the DNA fragments from 

1^^^ the gel matrix are also described in Sambrook et al. 
□ Any other non-destructive method of isolating DNA fragments of the desired size 

1 2 niay be employed. For example, size-based chromotography, HPLC, dHPLC or a sucrose 

density gradient could be used to reduce the DNA pool to those fragments within a 
1^ particular size range and then this smaller pool could be run on an electrophoresis gel. 

20 After isolation, adaptor sequences are ligated to the fragments. (Figure 2, Step 4) 

Adaptor sequences are generally oligonucleotides of at least 5 or 10 bases and preferably 
no more than 50 or 60 bases in length, however, adaptor sequences may be even longer, up 
to 100 or 200 bases depending upon the desired result. For example, if the desked 
outcome is to prevent amplification of a particular fragment, longer adaptor sequences 

25 designed to form stem loops or other tertiary structures may be ligated to the fragment. 
Adaptor sequences may be synthesized using any methods known to those of skill in the 
art. For the pxirposes of this invention they may, as options, comprise templates for PGR 
primers and/or tag or recognition sequences. The design and use of tag sequences is 
described in US Patent No. 5,800,992 and US Provisional Patent Application No. 


60/140,359, filed 6/23/99. Both of which are incorporated by reference for all purposes. 
Adaptor sequences may be Ugated to either blunt end or sticky end DNA. Methods of 
ligation will be known to those of skill in the art and are described, for example, in 
Sambrook et al. Methods include DNase digestion to "nick" the DNA, ligation with 

5 ddNTP and the use of polymerase I to fill in gaps or any other methods described in the art. 

Further complexity reduction is achieved by adding a specific nucleotide on the 5' 
end of the PGR primer as illustrated in Figure 3. The specific nucleotide further reduces 
the complexity of the resulting DNA pool because only those fragments which have been 
isolated after restriction enzyme digestion and contain the complement of the specific 

10 nucleotide(s) incorporated in the PGR primer will be amplified. Figure 3 A depicts the 

results of hybridization to an array after enzyme digestion, ligation to an adaptor and PGR 

^0 amplification. Figs. 3B and 3G depict the results of hybridization to an array after enzyme 

g digestion, ligation to an adaptor and PGR amplification where the PGR primers 

incorporated specific nucleotides in the 5' end of the primer. In Fig. 3B the 5' and 3* 

IB primers have different specific nucleotides incorporated. In Fig. 3A the 5' and 3* primers 
have the same nucleotides incorporated. The level of complexity in the isolated pool can 

;J be varied depending upon the identity and number of nucleotides incorporated into the 

m PGR primers. A number of embodiments of the present invention involve amplification by 

H PGR. Any of these embodiments may be further modified to reduce complexity using the 

20 above disclosed technique. 

Various methods of conducting PGR amplification and primer design and 
construction for PGR amplification will be known to those of skill in the art. PGR is a 
method by which a specific polynucleotide sequence can be amplified in vitro. PGR is an 
extremely powerful technique for amplifying specific polynucleotide sequences, including 

25 genomic DNA, single-stranded cDNA, and mRNA among others. As described in U.S. 
Pat. Nos. 4,683,202, 4,683,195, and 4,800,159 (which are incorporated herein by 
reference), PGR typically comprises treating separate complementary strands of a target 
nucleic acid with two oligonucleotide primers to form complementary primer extension 
products on both strands that act as templates for synthesizing copies of the desired nucleic 
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acid sequences. By repeating the separation and synthesis steps in an automated system, 
essentially exponential duplication of the target sequences can be achieved. Standard 
protocols may be found in^ for example Sambrook et al. which is hereby incorporated by 
reference for all purposes. 

5 In another embodiment, schematically illustrated in Figure 4, the step of 

complexity management of the DNA samples comprises digestion with a Type lis 
endonuclease thereby creating sticky ends comprised of random nucleic acid sequences, 
(Fig 4, Step 1) Type-IIs endonucleases are generally commercially available and are well 
known in the art. A description of Type lis endonucleases can be found in US Patent No. 
10 5,710,000 which is hereby incorporated by reference for all purposes. Like their Type-II 

, S counterparts, Type-IIs endonucleases recognize specific sequences of nucleotide base pairs 
within a double stranded polynucleotide sequence. Upon recognizing that sequence, the 

-p endonuclease will cleave the polynucleotide sequence, generally leaving an overhang of 

ui one strand of the sequence, or "sticky end." 

1^0 Type-II endonucleases, however, generally require that the specific recognition site 

P be palindromic. That is, reading in the 5' to 3' direction, the base pair sequence is the same 
, T for both strands of the recognition site. For example, the sequence 
m G4-A-A-T-T-C 
P C-T-T-A-A4-G 

20 is the recognition site for the Type-II endonuclease EcoRI, where the arrows indicate the 

cleavage sites in each strand. This sequence is palindromic in that both strands of the 

sequence, when read in the 5' to 3' direction are the same. 

The Type-IIs endonucleases, on the other hand, generally do not require 

palindromic recognition sequences. Additionally, these Type-IIs endonucleases also 
25 generally cleave outside of their recognition sites. For example, the Type-IIs endonuclease 

Earl recognizes and cleaves in the foUovnng manner: 

CTCTTCNtNNNN 

GAG AAG n n n n ^ n 
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where the recognition sequence is -C-T-C-T-T-C-, N and n represent complementary, 
ambiguous base pairs and the arrows indicate the cleavage sites in each strand. As the 
example illustrates, the recognition sequence is non-palindromic, and the cleavage occurs 
outside of that recognition site. 
5 Specific Type-IIs endonucleases which are useful in the present invention include, 

e.g., Earl, Mnll, Plel, Alwl, Bbsl, Bsal, BsmAI, BspMI, Esp3I, Hgal, Sapl, SfaNI, Bbvl, 
BsmFI, Fokl, BseRI, HphI and MboIL The activity of these Type-IIs endonucleases is 
illustrated in FIG. 5, which shows the cleavage and recognition pattems of the Type-IIs 
endonucleases. 

10 The sticky ends resulting from Type-IIs endonuclease digestion are then ligated to 

adaptor sequences (Fig 4, Step 2) Those of skill in the art v^U be familiar vsdth methods of 
^2 ligation. Standard protocols can be found in, for example, Sambrook et al., hereby 
M incorporated by reference for all purposes. Only those fragments containing the adaptor 

J 7^ sequence are isolated. (Figure 6) 

f P In addition to those methods of isolation discussed above, methods of isolation 

n which take advantage of unique tag sequences which may be constructed in the adaptor 

sequences may be employed. These tag sequences may or may not be used as PGR primer 
n templates. Fragments containing these tags can then be segregated from other non-tag 
i2 bearing sequences using various methods of hybridization or any of the methods described 
20 in the above referenced application. 

In another embodiment, depicted in Figure 18, the method of complexity reduction 
comprises digesting the DNA sample with two different restriction enzymes. The first 
restriction enzyme is a frequent base cutter, such as MSE I which has a four base 
recognition site. The second restriction enzyme is a rare base cutter, such as Eco RI, which 
25 has a 6 base recognition site. This results in three possible categories of fragments; (most 
common) those which have been cut on both ends with the frequent base cutter, (least 
common) those which have been cut on both ends v^th the rare base cutter, and those 
which have been cut on one end with the frequent base cutter and on one end with the rare 
base cutter. Adaptors are ligated to the fragments and PGR primers are designed such that 
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only those fragments which fall into the desired category or categories are amplified. This 
technique, employed with a six base cutter and a four base cutter can reduce complexity 8- 
fold when only those fragments from the latter category are amplified. Other combinations 
of restriction enzymes may be employed to achieve the desired level of complexity. 

5 In another embodiment, the step of complexity management comprises removing 

repetitive sequences. Figure 10 depicts a schematic representation of this embodiment. 
The nucleic acid sample is first fragmented. (Figure 10, Step 1) Various methods of 
fragmenting DNA will be known to those of skill in the art. These methods may be, for 
example, either chemical or physical in nature. Chemical fragmentation may include 

10 partial degradation with a DNAse, partial depurination with acid, the use of restriction 

enzymes or other enzymes which cleave DNA at known or unknown locations. Physical 
fragmentation methods may involve subjecting the DNA to a high shear rate. High shear 
rates may be produced, for example, by moving DNA through a chamber or channel with 

^ =f pits or spikes, or forcing the DNA sample through a restricted size flow passage, e.g., an 
aperture having a cross sectional dimension in the micron or submicron scale. 

In a preferred embodiment adaptor sequences are ligated to the resulting fragments. 
(Figure 10, Step 2) The fragments with or without adaptor sequences are then denatured. 

m (Figure 10, Step 3) Methods of denaturation will be will known to those of skill in the art. 

J^' After denaturation, the fragments are then allowed to reanneal. (Figure 10, Step 4) 

20 Annealing conditions may be altered as appropriate to obtain the level of repetitive 

sequence removal desired. Finally, double stranded sequences are removed (Figure 10, 
Step 4). Methods of removing double stranded sequences will be known to those of skill 
in the art and may include without limitation, methods of digesting double stranded DNA 
such as double strand specific nucleases and exonucleases or methods of physical 

25 separation including, without limitation gel based electrophoresis or size chromotography. 

In another embodiment, the step of complexity management comprises performing 
an arbitrarily primed polymerase chain reaction (AP PGR) upon the sample. AP PGR is 
described in US Patent No. 5,487,985 which is hereby incorporated by reference in its 
entirety for all purposes. Figure 7 depicts a schematic illustration of this embodiment. 



Performing AP PGR with random primers which have specific nucleotides incorporated 
into the primers produces a reduced representation of genomic DNA in a reproducible 
manner. Figure 8 shows the level of complexity reduction of human genomic DNA 
resulting from AP PGR with various primers. Golumn 1 lists the primer name. Golumn 2 

5 lists the primer sequence. Golumn 3 lists the annealing temperature. Golumn 4 lists the 

polymerase used. Golumn 5 lists the number correlated to a specific gene on the Hum6.8K 
GeneChip(R) probe array (Affymetrix, Inc. Santa Glara, Ga). Golumn 6 Usts the 
percentage of the human genes on the Hum6.8K GeneGhip(R) probe array found by 
fragments whose complexity has been reduced by this method. Figure 9 shows the 
10_ reproducibility of AP PGR. Independently prepared samples preps were subjected to AP 

i Q PGR using the same primers. The gel bands show that the level of reproducibility between 

Ji{ the samples is very high. 

P Primers may be designed using standard techniques. For example, a computer 

^ J program is available on the internet at the Operon Technologies, Inc. website at http: 
15^ www.operon.com. The Operon Oligo Toolkit allows a user to input a potential primer 
O sequence into the webfbrm. The site will instantly calculate a variety of attributes for the 
oligonucleotide including molecular weight, GG content, Tm, and primer-dimer sets. You 
may also plot the oligonucletoide against a second sequence. PGR amplification 
U techniques are described above in this application and will be well known to those of skill 
20 in the art. 

In another embodiment of the invention, the method reducing the complexity of a 
nucleic acid sample comprises hybridizing the sample to a nucleic acid probe containing a 
desired sequence which is bound to a solid support, such as a magnetic bead. For a 
description of hybridization of nucleic acids to solid supports, see US Pat No. 5,800,992 
25 incorporated by reference above. This sequence may comprise, for example, a sequence 
containing a SNP, a cDNA fragment, a chromosome fragment, a subset of genomic DNA 
or a subset of a library. The sequence may comprise as few as 16 nucleotides and may 
comprise as many as 2,000, 3,000, 5,000 or more nucleotides in length. Methods of 
designing and making oligonucleotide probes will be well known to those of skill in the 
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art. In one embodiment, the probe may contain a template sequence for a PGR primer. 
Solid supports suitable for the attachment of nucleic acid probe sequences will be well 
known to those of skill in the art but may include, glass beads, magnetic beads, and/or 
planar surfaces. Magnetic beads are commercially available jfrom, for example, Dynal 
5 (Oslo, Norway). The nucleic acid probes may be synthesized directly on the solid support 

or attached to the support as a full length sequence. Protocols for attaching magnetic beads 
to probes are included in US Patent No. 5,512,439 which is hereby incorporated by 
reference for all purposes. Standard hybridization protocols as discussed above may be 
employed. 

10 Figure 1 1 depicts a schematic representation of one example of the above 

,S embodiment, wherein the complexity management step is utilized to facilitate genome 
^ y Avide genotyping. Much of the cost of genotyping comes from multiplex PGR. In this 
.p embodiment, the entire sample preparation can be performed in a single tube without the 

need for multiplex PGR. Because the desired result is to genotype a DNA sample, the 
rfl desired sequence in Figure 1 1 contains a polymorphism. The oligonucleotide comprises 
O 32 bases with the SNP m the center. A magnetic bead is attached to the oligonucleotide 
3": probe. (Fig. 11, step 1) The probe is then exposed to, for example, fractionated genomic 
!U I>NA. (Fig.l 1, step 2). Adaptor sequences are ligated to both ends of the fragments. (Fig. 

1 1 11, step 3). The fragments are then amplified (Fig. 11, step 4) and the PGR product 

20 containing the desired polymorphism may then be analyzed by various methods including, 
for example, hybridization to an array or single base extension (SBE). SEE is described 
in, for example US Provisional Application 60/140,359 which is hereby incorporated by 
reference in its entirety for all purposes. 

The method may further comprise exposing the hybridized sample to a single 

25 strand DNA nuclease to remove the single stranded DNA. This embodiment may further 

comprise ligating an adaptor sequence containing a Glass II S restriction enzyme site to 
resulting duplexed DNA and digesting the duplex with the appropriate Glass II S 
restriction enzyme to release the attached sequences. The sequences are then isolated and 
a second adaptor sequence is ligated to the complex and the sequences are amplified. 
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Figures 12 and 13 depict schematic representations of an embodiment comprising 
the use of ClassIIs endonucleases. Both figures depict methods which may be employed 
for single tube genotyping without the need for multiplex PGR. In Figures 12 and 13, the 
desired sequence is a SNP. The oligonucleotide probe in Figure 12 is 32 bases long and in 

5 Figure 13 is 17 bases long. In both figures the SNP is in the center of the oligonucleotide. 
The oligonucleotide probe is bound to a magnetic bead. (Figs. 12 and 13, step 1). The 
probe is then hybridized to fragmented genomic DNA (Figs. 12 and 13, step 2). Single 
stranded DNA is digested with a single strand DNA nuclease leaving a DNA duplex 
attached to the magnetic bead, (Figs. 12 and 13, step 3). An adaptor sequence is then 

10 ligated to the duplex. The adaptor sequence contains a Class IIS restriction site. The 

probe length and Class IIS endonuclease are chosen such that the site where the duplex is 
:0 cut is between the SNP and the magnetic bead. In Figure 12 the Class IIS endonuclease 
S cuts directly adjacent to the SNP site, such that the SNP is part of the sticky end left by the 

H endonuclease digestion. (Fig. 12, step 5) In Figure 13 the endonuclease cuts closer to the 

1C= magnetic bead, leaving a number of bases between the sticky end and the SNP site. (Fig. 

L 13, step 5) In either case, the magnetic bead is released and the sequences are isolated. 

' 4 Ad^tor sequences are then ligated to the sticky ends. (Figs. 12 and 13, step 6) In both 

Figures 12 and 13 the adaptor sequences contain templates for PCR probes. The fragments 

5 "^ containing the SNP are then amplified (Figs. 12 and 13, step 7) and the PCR products may 

20 be analyzed in a nimiber of different methods including hybridization to an array designed 

to detect SNPs or SBE. 

In this embodiment, the adaptor sequence may further comprise a SNP 
identification sequence or tag. In this case, the array to which the PCR products are 
hybridized may be a generic tag array as described in the above referenced US Patent No. 

25 5,800,992 and US Provisional Patent Application 60/140,359 or a chimeric probe array 

(Figure 14). A chimeric probe array contains probes which interrogate both for particxxlar 
sequences characteristic of a genotype as well as for artificial sequences which have been 
ligated to specific fragments in the sample pool. This allows for higher specificity of 
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hybridization and better differentiation between probes. This embodiment is depicted in 
Figure 15. 

In another embodiment, depicted in Figure 16 the method of complexity reduction 
comprises hybridizing the DNA sample to a mismatch binding protein. Fig. 16, step 2. 
5 Mismatch binding proteins are described in Wagner R. and Radman, M. (1995) "Methods: 
A Companion to Methods in Enzymology" 7, 199-203 which is hereby incorporated by 
reference in its entirety for all purposes. Mismatch binding proteins preferentially bind to 
DNA duplexes which contain sequence mismatches. This allows for a relatively simple 
and rapid method to locate and identify SNPs. In this embodiment no prior knowledge of 

10 the SNP is required. Mismatch binding proteins are commercially available through 
GeneCheck (Ft. Collins, Co.). In a further embodiment, depicted in Figure 17, magnetic 

^0 beads are attached to the mismatch binding proteins. Mismatch binding proteins attached 
g to magnetic beads are commercially available through GeneCheck (Ft. Collins, Co.). After 
^7: hybridization the sample is digested with a 3' to 5' exonuclease (Fig. 16, step 3). 

11 Remaining single stranded DNA is then removed with a nuclease (Fig. 1 6, step 4). 

If it is desired to cut the duplex at the mismatch, then the enzyme resolvase may be 
I ";^ used. See US Patent Nos. 5,958,692, 5,871,91 1 and 5,876,941 (each of which is 
m incorporated by reference in their entireties for all purposes) for a description of varioxxs 
J3: methods of cleaving nucleic acids. The resolvases (e.g. X-solvases of yeast and 
20 bacteriophage T4, Jensch et al. EMBO J. 8, 4325 (1989)) are nucleolytic enzymes capable 
of catalyzing the resolution of branched DNA intermediates (e.g., DNA cruciforms) which 
can involve hundreds of nucleotides. In general, these enzymes are active close to the site 
of DNA distortion (Bhattacharyya et al, J. MoL Biol, 221, 1 191, (1991)). T4 
Endonuclease VII, the product of gene 49 of bacteriophage T4 (Kleff et al. The EMBO J. 
25 7, 1527, (1988)) is a resolvase (West, Annu. Rev. Biochem. 61, 603, (1992)) which was 

first shown to resolve Holliday-stractures (Mizuuchi et al.. Cell 29, 357, (1982)). T4 
Endonuclease VII has been shown to recognize DNA cruciforms (Bhattacharyya et al., 
supra; Mizuuchi et al., supra) and DNA loops (Kleff et al., supra), and it may be involved 
in patch repair. Bacteriophage T7 Endonuclease I has also been shown to recognize and 
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cleave DNA cruciforms (West, Ann. Rev. Biochem. 61, 603, (1992)). Eukaryotic 
resol vases, particularly from the yeast Saccharomyces cerevisiae, have been shown to 
recognize and cleave cruciform DNA (West, supra; Jensch, et al., EMBO J. 8, 4325 
(1989)). Other nucleases are known which recognize and cleave DNA mismatches. For 
5 example, SI nuclease is capable of recognizing and cleaving DNA mismatches formed 

when a test DNA and a control DNA are annealed to form a heteroduplex (Shenk et al., 
Proc. Natl. Acad. Sci. 72, 989, (1975)). The Nut Y repair protein of E. coli is also capable 
of detecting and cleaving DNA mismatches. 

10 Computer Implemented Analysis 

■f In another embodiment a computer system is used to model the reactions discussed 

^0 above to aid the user in selecting the correct experimental conditions. In this embodiment, 
^1 the sequence of the DNA sample must be known. A computer program queries an 
H electronic database containing the sequence of the DNA sample looking for sites which 
1# will be recognized by the enzyme being used. The method of modeling experiments can 

be employed for a wide variety of experiments. 
^ In one embodiment, the user can run multiple experiments altering various 

lU conditions. For example, if a user desires to isolate a particular sequence of interest in a 
ri fragment which has been digested with a restriction enzyme, the user can have the 
20 computer model the possible outcomes using a wide variety of restriction enzymes. The 

particular sequence which is selected may be chosen by specific criteria, i.e. because the 
region is beUeved to be associated with specific genes, polymorphisms, or phenotypes for 
example, or may be chosen at random. The user can then select the restriction enzyme 
which, for example, isolates the desired sequence in a fragment of unique size. 
25 Additionally or altematively, if the user desires to reduce complexity using the type IIS 

nuclease/ligation technique described above, the user can experiment v^th the length and 
sequence of the adaptors to determine the optimal sequence for the adaptors' "sticky" ends. 
This enables the user to be confident that they will obtain a fragment containing a 
particular sequence of interest or to fine tune the level of complexity in the DNA pool In 
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another embodiment, a user could model the kinetics of the denaturing, reannealing 
technique for removal of repeated sequences discussed above to determine the conditions 
Mrhich allow for the desired result. For example, a user may desire the removal of only a 
certain percentage of repeated sequences, 
5 For example, virtual restriction digests may be performed by querying an electronic 

database which contains the sequence of DNA of interest. Because the database contains 
the nucleic acid sequence and restriction enzymes cut at known locations based on the 
DNA sequence, one can easily predict the sequence and size of fragments which will result 
from a restriction digest of the DNA. Ideally, restriction enzymes which produce no two 
10 fragments of the same or very similar size are desired. Combinations of restriction 
^ enzymes may be employed. Those of skill in the art will be familiar with electronic 
;0 databases of DNA sequences. GenBank, for example, contains approximately 

2,570,000,000 nucleic acid bases in 3,525,000 sequence records as of April 1999. A 
H computer program searches the electronic database for a sequence which suits the 
ISO requirements of the particular restriction enzyme. For example, the restriction enzyme 
= Eco RI recognizes the sequence GAATTC and will cut a DNA molecule between the G 

. ^ and the first A. The computer program will query the chosen sequence for any occurences 
FU of the sequence GAATTC and mark the site where the restriction enzyme will cut. The 
^ program v^U then provide the user with a display of the resulting fragments. 
20 Exhibit 1 is an example of a program to conduct this type of virtual enzyme 

digestion. Exhibit 2 is an example of a program to virtually model the ligation of two 
sequences to each other. 

In another embodiment, the method of modeling experiments in a computer system 
can be used to design probe arrays. A database may be interrogated for any desired 
25 sequence, for example, a polymorphism. Computer modeled reactions are then performed 
to help determine the method for isolating a fragment of DNA containing the sequence of 
interest. These methods may comprise any of the methods described above, alone or in 
combination. Arrays are then constructed which are designed to interrogate the resulting 
fragments. It is important to note that for the purpose of designing arrays, the virtual 
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reactions need not be performed flawlessly, since the arrays may contain hundreds of 
thousands of sequences. 

One embodiment of the invention relies on the use of virtual reactions to 
predetermine the sequence of chosen DNA fragments which have subjected to various 
5 procedures. The sequence information for the chosen fragments is then used to design the 

probes which are to be attached to DNA arrays. Arrays may be designed and manufactured 
in any number of ways. For example, DNA arrays may be synthesized directly onto a 
solid support using methods described in, for example US Patent Nos. 5,837,832, 
5,744,305 and 5,800,992 and W095/11995 herein incorporated by reference for all 
10 purposes. See also, Fodor et al, Science, 251:767-777 (1991), Pirrung et al, U.S. Pat. No. 

5,143,854 (see also PCX Application No. WO 90/15070) and Fodor et al., PCT Publication 
j| No. WO 92/10092 and U.S. Pat. No, 5,424,186, each of which is hereby incorporated in its 
.p entirety by reference for all purposes. Techniques for the synthesis of these arrays using 
Jt; mechanical synthesis methods are described in, e.g., U.S. Pat. No. 5,384,261, incorporated 
vfi herein by reference in its entirety for all purposes Briefly, 5,837,832 describes a tiling 
p method for array fabrication whereby probes are synthesized on a solid support. These 
, 2 arrays comprise a set of oligonucleotide probes such that, for each base in a specific 
lU reference sequence, the set includes a probe (called the "wild-type" or "WT" probe) that is 
(1 exactly complementary to a section of the sequence of the chosen fragment including the 
20 base of interest and four additional probes (called "substitution probes"), which are 

identical to the WT probe except that the base of interest has been replaced by one of a 
predetermined set (typically 4) of nucleotides. Probes may be synthesized to query each 
base in the sequence of the chosen fragment. Target nucleic acid sequences which 
hybridize to a probe on the array which contain a substitution probe indicate the presence 
25 of a single nucleotide polymorphism. Other applications describing methods of designing 
tiling arrays include: US Patent Nos. 5,858,659, and 5,861,242 each of which is 
incorporated by reference in its entirety for all purposes. In a similar manner, arrays could 
be constructed to test for a variety of sequence variations including deletions, repeats or 
base changes greater than one nucleotide, US Patent Nos. 5,593,839 and 5,856,101 (each 
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of which is incorporated by reference for all purposes) describe methods of using 
computers to design arrays and lithographic masks. 

The label used to detect the target sequences will be determined, in part, by the 
detection methods being applied. Thus, the labeling method and label used are selected in 
5 combination with the actual detecting systems being used. Once a particular label has been 

selected, appropriate labeling protocols will be applied, as described below for specific 
embodiments. Standard labeling protocols for nucleic acids are described, e.g., in Maniatis; 
Kambara, H, et al. (1988) BioTechnology 6:816-821; Smith, L. et al. (1985) Nuc. Acids 
Res. 13:2399-2412; for polypeptides, see, e.g., Allen G. (1989) Sequencing of Proteins and 
10 Peptides, Elsevier, N.Y., especially chapter 5, and Greenstein and Winitz (1961) Chemistry 
; J of the Amino Acids, Wiley and Sons, N.Y. Carbohydrate labeling is described, e.g., in 
^0 Chaplin and Kennedy (1986) Carbohydrate Analysis: A Practical Approach, IRL Press, 

Oxford. Other techniques such as TdT end labeling may likewise be employed. 
y% Techniques for labeling protocols for use with SBE are described in, e.g. US Provisional 
M Patent Application 60/140,359 which is incorporated by reference above. 
PI . Generally, when using a DNA array a quickly and easily detectable signal is 
Z";^ preferred. Fluorescent tagging of the target sequence is often preferred, but other smtable 
m labels include heavy metal labels, magnetic probes, chromogenic labels (e.g., 
J ^ phosphorescent labels, dyes, and fluorophores) spectroscopic labels, enzyme linked labels, 
20 radioactive labels, and labeled binding proteins. Additional labels are described in U.S. 

Pat Nos. 5,800,992 and 4,366,241, and published PCT Application WO 99/13319 which 
are incorporated herein by reference. 

The hybridization conditions between probe and target should be selected such that 
the specific recognition interaction, i.e., hybridization, of the two molecules is both 
25 sufficiently specific and sufficiently stable. See, e.g., Hames and Higgins (1985) Nucleic 

Acid Hybridisation: A Practical Approach, IRL Press, Oxford. These conditions will be 
dependent both on the specific sequence and often on the guanine and cytosine (GC) 
content of the complementary hybrid strands. The conditions may often be selected to be 
universally equally stable independent of the specific sequences involved. This typically 
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will make use of a reagent such as an alkylammonium buffer. See, Wood et al. (1985) 
"Base Composition-independent Hybridization in Tetramethylammonium Chloride: A 
Method for Oligonucleotide Screening of Highly Complex Gene Libraries," Proc. Natl. 
Acad. Sci. USA, 82:1585-1588; and Krupov et al. (1989) "An Oligonucleotide 
5 Hybridization Approach to DNA Sequencing," FEES Letters, 256:1 18-122; each of which 

is hereby incorporated herein by reference. An alkylammonium buffer tends to minimize 
differences in hybridization rate and stability due to GC content. By virtue of the fact that 
sequences then hybridize with approximately equal affinity and stability, there is relatively 
little bias in strength or kinetics of binding for particular sequences. Temperature and salt 
10 conditions along with other buffer parameters should be selected such that the kinetics of 
'Iz renaturation should be essentially independent of the specific target subsequence or 

oligonucleotide probe involved. In order to ensure this, the hybridization reactions will 
M usually be performed in a single incubation of all the substrate matrices together exposed 
n to the identical same target probe solution under the same conditions. The hybridization 
l€l conditions will usually be selected to be sufficiently specific such that the fidelity of base 
Q matching will be properly discriminated. Of course, control hybridizations should be 

included to determine the stringency and kinetics of hybridization. See for example, US 
jU Patent No. 5,871,928 which is hereby incorporated in its entirety for all purposes. 

Another factor that can be adjusted to increase the ability of targets to hybridize to 
20 probes is the use of nucleic acid analogs of PNAs in the probes. They can be built into the 

probes to create a more uniform set of hybridization conditions across the entire array. See 
US Patent Application No. 08/630,427 incorporated by reference above. 

The detection methods used to determine where hybridization has taken place will 
typically depend upon the label selected. Thus, for a fluorescent label a fluorescent 
25 detection apparatus will typically be used. Pirrung et al. (1992) U.S. Pat. No. 5,143,854 

and Ser. No. 07/624,120, now abandoned, (both of which are hereby incorporated by 
reference for all purposes) describe apparatus and mechanisms for scaiming a substrate 
matrix using fluorescence detection, but a similar apparatus is adaptable for other optically 
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detectable labels. See also, US Patent Nos. 5,578,832, 5,834,758, and 5,837,832 each of 
which is incorporated by reference in its entirety for all purposes. 

A variety of methods can be used to enhance detection of labeled targets bound to a 
probe attached to a solid support. In one embodiment, the protein MutS (from E. coli) or 

5 equivalent proteins such as yeast MSHl, MSH2, and MSH3; mouse Rep-3, and 

Streptococcus Hex- A, is used in conjunction with target hybridization to detect probe- 
target complex that contain mismatched base pairs. The protein, labeled directly or 
indirectly, can be added during or after hybridization of target nucleic acid, and 
differentially binds to homo- and heteroduplex nucleic acid. A wide variety of dyes and 
10_ other labels can be used for similar purposes. For instance, the dye YOYO-1 is known to 
bind preferentially to nucleic acids containing sequences comprising runs of 3 or more G 
residues. Signal amplification methods as described in US Patent Application No. 

= P 09/276,774 may likewise be used. 

Q Various methods of hybridization detection will be known to those of skill in the 

iP art. See for example, US Patent Nos. 5,578,832, 5,631,734, 5,744,305 and 5,800,992 each 
P of which is hereby incorporated in its entirety fox all purposes. 

Examples 

20 

Example 1 - Restriction Enzyme Digest/Sizing 

The complexity of total genomic DNA from human and yeast was reproducibly 
reduced using a restriction enzyme digestion. For each species 0.5 ug genomic DNA was 
digested with 20 units of EcoRI in a total volume of 40 ul at 37 °C overnight (Figure 2, 
25 Stepl). The enzyme was inactivated by incubation at 65 for 10 minutes. 

The DNA solution was mixed with 10 ul 5x loading buffer and separated by gel 
electrophoresis on a 2% agarose gel. (Figure 2, Step 2) The gel was visualized by 
ethidium bromide staining. Fragments of 250 - 350 bp were excised from the gel and 
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purified using a QIAquick gel extraction kit (Qiagen). (Figure 2^ Step 3) Alternatively, 
fragments of the required size could have been isolated using HPLC. 

Adaptor sequences containing PGR primer template sequences were then ligated to 
the purified fragments using lOOU T4 ligase in Ix T4 DNA ligase buffer (New England 
5 Biolabs) at 16 overnight. The adaptor sequences were 5'- 

d(pAATTCGAACCCCTTCGGATC)-3' and 5'-d(GATCCGAAGGGGTTCGAATT)-3' 
(Figure 2, Step 4) The ligase was then heat inactivated at 65 for 15 minutes. 

The fragments were then subjected to PGR with one primer that corresponded to 
the PGR primer template sequence 5'-d(GATGGGAAGGGGTTCGAATT)-3' (Figure 2, 
10,^, Step 5). The PGR mixture contained approx. 1 ng ligated DNA fragments, 5 units 

AmpliTaq Gold polymerase (Perkins Elmer), 5 uM pimer, 200uM dNTPs, 15 mM Tris- 
Q HCl (pH8.2), 50 mM KGl, 2,5 mM MgGl2 in a final volume of 50 ul. PGR was performed 
;£ in a Perkin-Elmer 9600 thermocycler using an initial 10 minute denaturation at 95 °G, 35 
cycles of a 1 minute denaturation at 94 °G, annealing for 1 minute at 57 and extension 
15^ at 72 for 2 minutes. This is followed by a final 5 minute extension cycle at 72 °G. 

The PGR products were then purified with QIAquick PGR Purification kit (Qiagen) 
according to the manufacturer's instructions and fi-agmented with DNase L 
O The remaining fragments were then labeled with biotin-N6-ddATP as follows: In 

each tube, incubate 10 ug DNA with 0.3 unit Dnasel (Promega) at 37 "^G for 30 minutes in 
20 a 45 ul mixture also containing 10 mM Tris-Actate (pH 7.5), 10 mM magnesium acetate 
and 50 mM potassium acetate. Stop the reaction by heating the sample to 95 °G for 15 
minutes. Label the sample by adding 60 unit terminal transferase and 4 pmol biotin-N6- 
ddATP (Dupont NEN) followed by incubation at 37 *^G for 90 minutes and a final heat 
inactivation at 95 °C for 15 minutes. 
25 The labeled DNA was then hybridized to an array in a hybridization mixture 

containing 80 ug labeled DNA, 160 ug human GOT-1 DNA (GIBGO), 3.5 M 
tetramethylamonium cloride, 10 mM MES (pH 6.5), 0.01% Triton-100, 20 ug herring 
sperm DNA, 100 ug bovine serum albumin and 200 pM control oligomer at 44 °C for 40 
hours on a rotisserie at 40 rpm. The arrays were then washed with 0.1 M MaGl in 10 mM 



MES at 44 for 30 minutes on a rotisserie at 40 rpm. The hybridized arrays were then 
stained with a staining solution [10 mM MES (ph 6.5), 1 M NaCl, 10 ug/ml steptaviden R- 
phycoerythrin, 0.5 mg/ml acetylated BSA, 0.01% Triton-100] at 40 °C for 15 minutes. 
The arrays were then washed with 6x SSPET [0.9 M NaCl, 60 mM NaH2P04 (pH 7.4), 6 

5 mM EDTA, 0.005 % Triton- 100] on a GeneChip® Fluidics Station (Affymetrix, Inc., 

Santa Clara, CA) 10 times at 22 ""C. The arrays were then anti-streptavidin antibody stained 
at 40 ""C for 30 minutes with antibody solution [lOmM MES (pH 6.5), 1 M NaCl, 10 ug/ml 
streptavidin R-phycoerythrin, 0.5 mg/ml actylated BSA, 0,01% Triton-lOO], The arrays 
are then restained with staining solution for 15 minutes followed by 6X SSPET washing as 
ipi^ above. The arrays are then scanned with a confocal scanner at 560 nm. The hybridization 

^0 patterns were then screened for SNP detection with a computer program as described in 

□ D.G. Wang et al Science 280, 1077-1082, 1998. The results of the hybridization can be 

% seen in Figures 8 A and 8B. 


1^ Example 2 - Digestion with a Type lis Endonuclease and Selective Ligation 
J 5 Complexity was reproducibly reduced after digestion with a type lis endonuclease 

1^ and selective ligation to an adaptor sequence. 2 ug of genomic I>NA was digested with 
20 Bbv I at 37 ovemight. (Figure 3, Step 1) The enzyme was heat inactivated at 65 for 
15 minutes. 

Adaptors containing PGR primer template sequences were ligated in a 50 ul 
mixture of 400 ng digested genomic DNA, 10 pmol adaptor and 40 unit T4 ligase in a IX 
T4 ligase buffer. (Figure 3, Step 2) The adaptor sequences were as follows: 5'- 
25 d(pATNNGATCCGAAGGGTTCGAATTC)-3 ' and 

5'GAATTCGAACCCCTTCGGATC)-3'. The ligation was conducted at le^'C overnight. 
The ligase was inactivated by incubation at 65°C for 15 minutes. 

The fragments were then subjected to PGR with one primer that corresponded to 
the PGR primer template sequence: 5'-GAATTCGAACCCCTTCGGATC)-3' in a 50 ul 
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reaction containing 20 ng ligated DNA, 1 unit AmpliTaq Gold polymerase (Perkins 
Elmer), 3 uM primer, 200uM dNTPs, 15 mM Tris-HCl (pH8.0), 50 mM KCl, 2.5 mM 
MgCb. PGR was performed in a Perkin-Elmer 9600 thermocycler using an initial 10 
minute denaturation at 95^C, 35 cycles of a 0.5 minute denaturation at 94°C, aimealing for 
5 0.5 minute at 57°C and extension at 72^C for 2 minutes. This is followed by a final 5 
minute extension cycle at 72°C. 

Example 3 - Double Digestion and Selective PGR 

Human genomic DNA was digested in a 40 ul reaction at 37 °C for 1 hour. The 
lf>j reaction mixture contained 0.5 ug human genomic DNA, 0.5 mM DTT, 5 unit EcoRI (New 
^3 England Biolabs), 5 units Sau3AI (New England Biolabs), 0.5 ng/ul BSA, 10 mM Tris- 
O Acetate (pH 7.5), 10 mM magnesium acetate and 50 mM potassiimi acetate. The enzymes 
J 5 were inactivated at 65 °G for 15 minutes. 

The restriction fragments were then ligated to adaptor sequences. The ligation 
i f mixture contained: 5 pmol Eco R I adaptor [5'-d(pAATTCGAAGGCCTTCGGATC)-3' 
i:j and 5'-d(GATCGGAAGGGGTTCG)-3'], 50 pmol Sau3A I adaptor [5'- 

d(pGATCGCCCTATAGTGAGTCGTATTACAGTGGACCATCGAGGGTCA).3'],5 
O mM DTT, 0.5 ng/ul BSA, 100 unit T4 DNA ligase, 1 mM ATP, 10 mM Tris-Acetate (pH 
7.5), 10 mM magnesium acetate and 50 mM potassium acetate]. The ligation mixture was 
20 incubated v^th the restriction fragments at 37^C for 3 hours. The ligase was inactivated at 
65 ""C for 20 minutes. 

The ligated DNA target was then amplified by PGR. The PGR mixture contained 
12.5 ng ligated DNA, 1 imit AmpliTaq Gold polumerase (Perkins Elmer), 0.272 mM 
EcoRI selective primer (5'-AAGGGGTTGGGAATTCGG-3'; GG as the selective bases), 
25 0.272 uM Sau3AI selective primer (5'-TGAGTATAGGGGGATGTG-3'; TG as the 

selective bases), 200 uM dNTPs, 15 mM Tris-HCl (pH 8,0), 50 mM KCl, 2.5 mM MgCk 
in a final volume of 50 ul. PGR was performed in a Perkin-Elmer 9600 thermocycler using 
an initial 10 minute denaturation at 95 ^G, 35 cycles of a 1 minute denaturation at 94 °G, 


annealing for 1 minute at 56 ""C and extension at 72 for 2 minutes. This is followed by a 
final 5 minute extension at 72 °C, 

Example 4, Arbitarily Primed PGR 
5 PGR pimers were designed with the Operon Oligo Toolkit described in the 

specification above. 

Human genomic DNA was amplified in a 100 ul reaction containing 100 ng 

genomic DNA, 1 .25 units AmpliTaq Gold polymerase (Perkin Elmer), 10 uM arbitary 

primer, 200 mM dNTPs, 10 mM tris-HCl (pH 8.3), 50 mM KGl and 2.5 mM MgCh. 
10„^ PGR was performed in a Perkin-Elmer 9600 thermocycler using an initial 10 

^0 minute denaturation at 95 °G, 35 cycles of a 1 minute denaturation at 94 °G, annealing for 
i3 1 minute at 56 ^G and extension at 72 for 2 minutes. This is followed by a final 7 minute 
% extension at 72 ^G. 

W The PGR product was then purified, fi:agmented, labeled and hybridized as 

15 described in the examples above. 

Example 5 - SNP discovery - Generally 

As an example, the present invention may be directed to a method for simplifying 
? the detection of or comparing the presence of absence of SNPS among individuals, 

20 populations, species or between different species. This invention allows for a quick and 

cost-effective method of comparing polymorphism data between multiple individuals. 
First, a reduced representation of a nucleic acid sample is produced in a repeatable and 
highly reproducible maimer from multiple individuals, using any of the above described 
techniques alone or in combination. Then, the data generated by hybridizing the DNA 

25 samples collected fi-om multiple individuals to identical arrays in order to detect for the 

presence or absence of a number of sequence variants is compared. Arrays are designed to 
detect specific SNPS or simply to detect the presence of a region known to fi:equently 
contain SNPS, In the latter case, other techniques such as sequencing could be employed 
to identify the SNP. 
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SNP discovery - method 1 

Typically, the detection of SNPs has been made using at least one procedure in 
which the nucleic acid sequence that may contain the SNP is amplified using PGR primers. 

5 This use can create an expense if many SNPs are to be evaluated or tested and it adds 

significantly more time to the experiment for primer design and selection and testing. The 
following example eliminates the need for the specific PGR amplification step or steps. 
First, using the methods provided in example 1 above, a restriction enzyme or enzymes is 
used to cut genomic DNA at a large number of sites and a size range of restriction 
10 fragments is selected for assay. An electronic database, such as GenBank is queried to 

determine which sequences would be cut v^th the specific restriction enzyme(s) that were 

' 0 selected above. The sequences of the resulting fragments are then used to design DNA 

arrays which will screen the regions for the SNPs or other variants. The selected fragments 
are then subjected to further fragmentation and hybridized to the array for analysis. 

IfO 

SNP discovery - Method 2 
Z;^ Alternatively, the method provided in example 2 above may be employed, type IIS 

1 U restriction enzymes cut genomic DNA from each individual and adaptor sequences are 
l2 designed to ligate to specific fragments as desired. Adaptor sequences may include both 
20 random and specific nucleotide ends as required to produce the desired result. If desired, 

amplification primers may be designed to hybridize to the adaptor sequences, allowing for 
amplification of only the fragments of interest. An electronic database and computer 
modeling system may be used to aid in the selection of appropriate experimental 
conditions and to design the appropriate arrays. The fragments are then hybridized to the 
25 array for analysis. 


SNP discovery - Method 3 

As another altemative, MutS Protein were used to isolate DNA containing SNPS 
for analysis on an array. 3 ugs of DNA was fragmented with Eco R I (altematively a 



Dnase I could have been used.) At this point an equal amount of control DNA was added 
(this step is optional). 

0.5ug of the fragments were denatured at 95 °C for 10 minutes and gradually 
cooled to 65 °C over a 60 minute period. The fragments were then incubated at 65 °C for 
30 minutes and the temperature was ramped down to 25 °C over a 60 minute period. 1 .5 
ug MutS protein (Epicentre) was then added and allowed to incubate at room temperature 
for 1 5 minutes to allow for binding. (Figure 7, Step 1) 

The bound fragments were then digested with 20 xmits T7 polymerase (New 
England Biolabs) at 30 ""C for 30 minutes. (Figure 7, Step 2) The T7 polymerase was 
inactivated by incubation at 65 for 10 minutes. 

Single stranded DNA was trimmed with 100 units of nuclease SI (Boehringer- 
Mannheim) at 16 '^C for 15 minutes. (Figure 7, Step 3) The enzymes inactivated by 
adding 50 nmol EDTA and incubation at 65°C for 15 minutes. 

Adaptor sequences containing PGR primer templates were then ligated to the DNA 
sequences in a 10 ul ligation mixture: lul DNA solution, 4 ul dH20, 1 ul lOX T4 DNA 
ligase buffer, 3 ul 10 mM adaptor [5'-d(GATCCGAAGGGGTTCGAATT)-3' and 5'- 
d(pGAATTCGAACCCCTTCGGATC-eO and 1 ul 400 U/ul T4 DNA ligase] and 
incubated at 16 °C overnight and then inactivated at 65 °C for 15 minutes. (Figure 7, Step 
4) 

The sequences were amplified in a 25 ul reaction containing 0.25 pmol template 
DNA, 0.125 units AmpliTaq Gold polymerase (Perkin Elmer), 3 uM primer, [5'- 
d(GATCCGAAGGGGTTCGAATT)-3'], 200 uM dNTPs, 15 mM tris-HCl (pH 8.0), 50 
mM KCl and 1.5 mM MgCls. 

PGR was performed in a MJ Research Tetrad thermocycler using an initial 10 
minute denaturation at 95 °C, 35 cycles of a 0.5 minute denaturation at 94 °C, annealing 
for 0,5 minute at 57 and extension at 72 °C. This is followed by a final 5 minute 
extension at 72 *^C. 

The sequences were then labeled and hybridized to an array as described above. 



SNP discovery - Method 4 

As another alternative, oligonucletides attached to magnetic beads may be used for 
allele specific SNP enrichment and genotyping. Synthesized biotin-tagged 
oligonucleotides containing sequences complementary to the regions of desired SNPs 
were mixed with target DNA in a 1000: 1 ratio. (Alternatively, a 10:1, 20:1, 50:1, 250:1 or 
any other ratio could have been chosen.) 

The sample was then denatured at 95 °C for 10 minutes allowed to reanneal by 
slowly cooling to room temperature. 

The sample was then bound to streptavadin-magnetic beads (Promega) by mixing 
the sample and the beads and incubation at room temperature for 10 minutes. The beads 
were then washed with IX MES with IM Sodium Chloride (NaCl) three times. The beads 
were then resuspended in 50 ul IX mung bean nuclease buffer and mixed vdth 1 unit of 
mung bean nuclease. The beads were then incubated at 30*^C for 15 minutes. The mung 
bean nuclease was then inactivated by adding 1% SDS. The beads were then washed with 
IX MES with IM NaCl three times. 

The beads were then resuspended in ligation mixture containing T4 ligase in 1 X 
T4 ligase buffer and 200 fold excess adaptor I sequence [5'- 
d(ATTAACCCTCACTAAAGCTGGAG)-3'and 5'- 

d(pCTCCAGCTTTAGTGAGGGTTAAT)-3' Bpml recognition sites are highlighted in 
boldface] at 16 °C overnight. The ligase was then inactivated by incubation at 65 *^C for 10 
minutes. 

The beads were then washed with IX MES with IM NaCl three times and then 
resuspended in 50 ul IX Bpm I restriction buffer. BPM I was then added and the beads 
were incubated at 37 °C for 1 hr. The enzyme was inactivated by incubation at 65 °C for 
10 minutes and the supernatant solution with the sequences containing the desired SNPs 
was collected. 

A second set of adaptor sequences containing PCR template sequences [5'- 
d(pCTATAGTGAGTCGTATT-3') and (5'-AATACGACTCACTATAGNN-3')] and 



ligase were then added to the supernatant solution and incubated at 16 °C overnight. The 
ligase was then heat inactivated at 65 °C for 10 minutes. 

The samples were then amplified with PCR using T3 (S'- 
ATTAACCCTCACTAAAG-S') and T7 5'-d(TAATACGACTCACTATAGGG)-3' 
5 sequencing primers (Operon) in a 50 ml reaction containing 10^ copies of each target 

DNA, 1 unit AmpliTaq Gold polymerase (Perkin Elmer), 2 uM each primer, 200 uM 
dNTPs, 15 mM tris-HCl (pH 8.0), 50 mM KCl and 2.5 mM MgC^. 

PCR was performed in a MJ Research Tetrad Thermocycler using an initial 10 
minute denaturation at 95 ^C, 45 cycles of a 0.5 minute denaturation at 94 °C, aimealing 
10 for 0.5 minute at 52 °C and extension at 72 °C for 1 minute. This is followed by a final 5 

minute extension at 72 *^C. The fragments were then labeled and hybridized to an array. 

=P Methods of Use 

Id The present methods of sample preparation and analysis are appropriate for a wide 

15^"' variety of applications. Any analysis of genomic DNA may be benefitted by a 
\ 3 reproducible method of complexity management. 

As a preferred embodiment, the present procedure can be used for SNP discovery 
and to genotype individuals. For example, any of the procedures described above, alone 
I J. or in combination, could be used to isolate the SNPs present in one or more specific 
20 regions of genomic DNA. Arrays could then be designed and manufactured on a large 
scale basis to interrogate only those fragments containing the regions of interest. 
Thereafter, a sample from one or more individuals would be obtained and prepared using 
the same techniques which were used to design the array. Each sample can then be 
hybridized to a pre-designed array and the hybridization pattem can be analyzed to 
25 determine the genotype of each individual or a population of individuals as a whole. 
Methods of use for polymorphisms can be found in, for example, co-pending U.S. 
application 08/813,159. Some methods of use are briefly discussed below. 


Correlation of Polymorphisms with Phenotypic Traits 
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Some polymorphisms occur within a protein coding sequence and contribute to 
phenotype by affecting protein structure. The effect may be neutral, beneficial or 
detrimental, or both beneficial and detrimental, depending on the circumstances. For 
example, a heterozygous sickle cell mutation (which involves a single nucleotide 
polymorphism) confers resistance to malaria, but a homozygous sickle cell mutation is 
usually lethal. Other polymorphisms occur in noncoding regions but may exert phenotypic 
effects indirectly via influence on replication, transcription, and translation. A single 
polymorphism may affect more than one phenotypic trait. Likewise, a single phenotypic 
trait may be affected by polymorphisms in different genes. Further, some polymorphisms 
predispose an individual to a distinct mutation that is causally related to a certain 
phenotype. 

Phenotypic traits include diseases that have known but hitherto unmapped genetic 
components (e.g., agammaglobulimenia, diabetes insipidus, Lesch-Nyhan syndrome, 
muscular dystrophy, Wiskott-Aldrich syndrome, Fabry's disease, familial 
hypercholesterolemia, polycystic kidney disease, hereditary spherocytosis, von 
Wiliebrand's disease, tuberous sclerosis, hereditary hemorrhagic telangiectasia, familial 
colonic polyposis, Ehlers-Danlos syndrome, osteogenesis imperfecta, and acute 
intermittent porphyria). Phenotypic traits also include symptoms of, or susceptibility to, 
multifactorial diseases of which a component is or may be genetic, such as autoimmune 
diseases, inflammation, cancer, diseases of the nervous system, and infection by 
pathogenic microorganisms. Some examples of autoimmune diseases include rheumatoid 
arthritis, multiple sclerosis, diabetes (insulin-dependent and non-independent), systemic 
lupus erythematosus and Graves disease. Some examples of cancers include cancers of the 
bladder, brain, breast, colon, esophagus, kidney, leukemia, liver, lung, oral cavity, ovary, 
pancreas, prostate, skin, stomach and uterus. Phenotypic traits also include characteristics 
such as longevity, appearance (e.g., baldness, obesity), strength, speed, endurance, fertility, 
and susceptibility or receptivity to particular drugs or therapeutic treatments. 

Correlation is performed for a population of individuals who have been tested for 
the presence or absence of a phenotypic trait of interest and for polymorphic markers sets. 



To perform such analysis, the presence or absence of a set of polymorphisms (i.e. a 
polymorphic set) is determined for a set of the individuals, some of whom exhibit a 
particular trait, and some of which exhibit lack of the trait. The alleles of each 
polymorphism of the set are then reviewed to determine whether the presence or absence 
of a particular allele is associated with the trait of interest. Correlation can be performed 
by standard statistical methods such as a K-squared test and statistically significant 
correlations between polymorphic form(s) and phenotypic characteristics are noted. For 
example, it might be found that the presence of allele Al at polymorphism A correlates 
with heart disease. As a further example, it might be found that the combined presence of 
allele Al at polymorphism A and allele Bl at polymorphism B correlates with increased 
milk production of a farm animal. (See, Beitz et al., US 5,292,639 

Genetic Mapping of Phenotypic Traits 

Linkage analysis is useful for mapping a genetic locus associated with a phenotypic 
trait to a chromosomal position, and thereby cloning gene(s) responsible for the trait. See 
Lander et al., Proc. Natl. Acad. Sci. (USA) 83, 7353-7357 (1986); Lander et al., Proc. Natl. 
Acad. Sci. (USA) 84, 2363-2367 (1987); Donis-Keller et al.. Cell 51, 319-337 (1987); 
Lander et al.. Genetics 121, 185-199 (1989)). Genes localized by Imkage can be cloned by 
a process known as directional cloning. See Wainwright, Med. J. Australia 159, 170-174 
(1993); Collins, Nature Genetics 1, 3-6 (1992) (each of which is incorporated by reference 
in its entirety for all purposes). 

Lmkage studies are typically performed on members of a family. Available 
members of the family are characterized for the presence or absence of a phenotypic trait 
and for a set of polymorphic markers. The distribution of polymorphic markers in an 
mformative meiosis is then analyzed to determine which polymorphic markers co- 
segregate with a phenotypic trait. See, e.g., Kerem et al., Science 245, 1073-1080 (1989); 
Monaco et al., Nature 316, 842 (1985); Yamoka et al.. Neurology 40, 222-226 (1990); 
Rossiter et al., FASEB Journal 5, 21-27 (1991). 



Disequilibrium mapping of the entire genome 

Linkage disequilibrium or allelic association is the preferential association of a 
particular allele or genetic marker with a specific allele, or genetic marker at a nearby 
chromosomal location more frequently than expected by chance for any particular allele 
5 frequency in the population. For example, if locus X has alleles a and b, which occur 

equally frequently, and linked locus Y has alleles c and d, which occur equally frequently, 
one would expect the combination ac to occur with a frequency of 0,25, If ac occurs more 
frequently, then alleles a and c are in linkage disequilibrium. Linkage disequilibrium may 
result from natural selection of certain combination of alleles or because an allele has been 
10 introduced into a population too recently to have reached equilibrium with linked alleles, 
;S A marker in linkage disequilibrium can be particularly usefiil in detecting 

^2 susceptibility to disease (or other phenotype) notwithstanding that the marker does not 
i cause the disease. For example, a marker (X) that is not itself a causative element of a 

11 disease, but which is in linkage disequilibrium with a gene (including regulatory 

isO sequences) (Y) that is a causative element of a phenotype, can be detected to indicate 
Q susceptibility to the disease in circumstances in which the gene Y may not have been 
y identified or may not be readily detectable. 

12 Marker assisted breeding 

20 Genetic markers can decipher the genomes in animals and crop plants. Genetic 

markers can aid a breeder in the understanding, selecting and managing of the genetic 
complexity of an agronomic or desirable trait. The agriculture world, for example, has a 
great deal of incentive to try to produce food with a rising number of desirable traits (high 
yield, disease resistance, taste, smell, color, texture, etc.) as consumer demand and 

25 expectations increase. However, many traits, even when the molecular mechanisms are 
known, are too difficult or costly to monitor during production. Readibly detectable 
polymorphisms which are in close physical proximity to the desired genes can be used as a 
proxy to determine whether the desired trait is present or not in a particular organism. This 
provides for an efficient screening tool which can accelerate the selective breeding process. 



Pharmacogenomics 

Genetic information can provide a powerful tool for doctors to determine what 
course of medicine is best for a particular patient. A recent Science paper entitled 
"Molecular Classification of Cancer: Class Discovery and Class Prediction by Gene 
Expression Monitoring/' (to be published 10/15/99 hereby incorporated by reference in its 
entirety for all purposes) discusses the use of genetic information discovered through the 
use of arrays to determine the specific type of cancer a particular patient has. The paper 
goes on to discuss the ways in which particular treatment options can then be tailored for 
each patient's particular type of cancer. Similar uses of genetic information for treatment 
plans have been disclosed for patients with HIV. (See US Patent Application 5,861,242). 

The pharmaceutical industry is likevsdse interested in the area of 
pharmacogenomics. Every year pharmaceutical companies suffer large losses from drugs 
which fail clinical trials for one reason or another. Some of the most difficult are those 
drugs which, while being highly effective for a large percentage of the population, prove 
dangerous or even lethal for a very small percentage of the population. Pharmacogenomics 
can be used to correlate a specific genotype with specific responses to a drug. The basic 
idea is to get the right drug to the right patient. If pharmaceutical companies (and later, 
physicians) can accurately remove from the potential recipient pool those patients who 
would suffer adverse responses to a particular drug, many research efforts which are 
currently being dropped by pharmaceutical companies could be resurrected saving 
hundreds of thousands of dollars for the companies and providing many currently 
imavailable medications to patients. 

Similarly, some medications may be highly effective for only a very small 
percentage of the population while proving only slightly effective or even ineffective to a 
large percentage of patients. Pharmacogenomics allows pharamaceutical companies to 
predict which patients would be the ideal candidate for a particular drug, thereby 
dramatically reducing failure rates and providing greater incentive to companies to 
continue to conduct research into those drugs. 



Forensics 

The capacity to identify a distinguishing or unique set of forensic markers in an 
individual is useful for forensic analysis. For example, one can determine whether a blood 
sample from a suspect matches a blood or other tissue sample from a crime scene by 
determining whether the set of polymorphic forms occupying selected polymorphic sites is 
the same in the suspect and the sample. If the set of polymorphic markers does not match 
between a suspect and a sample, it can be concluded (barring experimental error) that the 
suspect was not the source of the sample. If the set of markers does match, one can 
conclude that the DNA from the suspect is consistent with that found at the crime scene. If 
frequencies of the polymorphic forms at the loci tested have been determined (e.g., by 
analysis of a suitable population of individuals), one can perform a statistical analysis to 
determine the probability that a match of suspect and crime scene sample would occur by 
chance. 

Patemity Testing / Determination of Relatedness 

The object of patemity testing is usually to determine whether a male is the father 
of a child. In most cases, the mother of the child is known and thus, the mother's 
contribution to the child's genotype can be traced. Patemity testing investigates whether 
the part of the child's genotype not attributable to the mother is consistent with that of the 
putative father. Patemity testing can be performed by analyzing sets of polymorphisms in 
the putative father and the child. Of course, the present invention can be expanded to the 
use of this procedure to determine if one individual is related to another. Even more 
broadly, the present invention can be employed to determine how related one individual is 
to another, for example, between races or species. 



Conclusion 

From the foregoing it can be seen that the advantage of the present invention is that 
it provides a flexible and scalable method for analyzing complex samples of DNA, such as 
genomic DNA. These methods are not limited to any particular type of nucleic acid 
sample: plant, bacterial, animal (including human) total genome DNA, RNA, cDNA and 
the like may be analyzed using some or all of the methods disclosed in this invention. This 
invention provides a powerful tool for analysis of complex nucleic acid samples. From 
experiment design to isolation of desired fragments and hybridization to an appropriate 
array, the above invention provides for faster, more efficient and less expensive methods of 
complex nucleic acid analysis. 

All publications and patent applications cited above are incorporated by reference 
in their entirety for all purposes to the same extent as if each individual publication or 
patent application Mrere specifically and individually indicated to be so incorporated by 
reference. Although the present invention has been described in some detail by v^ay of 
illustration and example for purposes of clarity and understanding, it Mdll be apparent that 
certain changes and modifications may be practiced within the scope of the appended 
claims. 


